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Abstract 

The simultaneous combination of scanning probe methods (tunnelling and force microscopies, 
STM and ATM) is a unique way to get an information about crystallographic and electronic 
structure of the studied surface. Here we apply these methods accompanied by the state-of-the- 
art density functional theory (DFT) calculations to shed a light on the structure and electronic 
properties of the strongly-corrugated graphene/Rh(lll) system. The atomically resolved images 
are obtained for both STM and AFM modes and compared with the DFT results showing good 
agreement between theory and experiment. 
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I. INTRODUCTION 


The discovery of the fascinating properties of graphene BE] initiates the fountain of 
experimental and theoretical works devoted to studies of its electronic and structural prop¬ 
erties. In most cases a graphene layer in these works is in contact with the supporting 
substrate: insulating, semiconducting or metallic. The latter supports are considered as the 
most interesting from technological and fundamental-science points of view and several re¬ 
cent reviews have pointed out the importance of these graphene/metal systems [3H6]. Firstly, 
the preparation of graphene on metals is presently considered as a cheapest way to produce 
large graphene layers of different thickness that can be transferred onto the polymer or insu¬ 
lating support and then used for preparation of different devices [M]. Another interesting, 
more fundamental, aspect of the graphene-metal interface is the nature of interaction at the 
interface |6]. Here, the intuitive description divides the whole set of the graphene-metal in¬ 
terfaces on two large subclasses of the weakly and strongly interacting graphene with metal. 
As a measure of interaction strength one can consider the bonding energy per carbon atom 
between graphene and metal, the doping level of graphene, or/and the state and position 
of the Dirac cone in the electronic structure of the system. Interesting to note that this 
“graphene-metal interaction” puzzle is valid for a graphene layer adsorbed on metal as well 
as for the metal layers deposited on graphene (free-standing or substrate-supported). 

Among the graphene/metal interfaces, the graphene layers on the lattice-mismatched 
close-packed surfaces [Rh(lll), Ru(OOOl), Ir(lll), Pt(lll)] are under intent attention [lUT 
[H]. The interaction between graphene and metal is abruptly changed from strong to weak 
interaction when going from Rh or Ru to Ir or Pt. Moreover, these systems are interesting 
from the perspective that metal atoms deposited on top of such graphene/metal systems will 
form an ordered arrays of clusters due to the moire structure of the graphene layer on the 
lattice-mismatched surfaces. In such moire structures the geometry of the system (buckling 
of a graphene layer) as well as the space variation of the charge (potential) will define the 
structure and the electronic properties of the metal-clusters/graphene/metal systems. 

Here we performed the studies of the crystallographic and electronic structure of the 
strongly-corrugated graphene layer on Rh(lll) via combination of scanning probe techniques 
(STM and AFM) and density-functional theory (DFT) calculations. These results show 
that a graphene layer in this system is strongly corrugated and electron density is strongly 
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localised at the graphene/Rh(lll) interface in places where graphene is strongly bonded 
to the substrate. Our results give also an information on the role of dispersive forces (van 
der Waals) in the bonding mechanism in this system and help to understand the clustering 
mechanism on top of strongly corrugated graphene layer. 

II. EXPERIMENTAL 

The electronic and structural properties of the graphene/Rh(lll) system are obtained 
using generalised gradient approximation as parameterized by Perdew et al. (PBE) [18] to 
the exchange correlation potential. For solving the resulting Kohn-Sham equation we use 
the Vienna Ab Initio Simulation Package (VASP) [T9| with the projector augmented wave 
basis sets [201 • The long-range van der Waals interactions were accounted for by means of 
a DFT-D2 approach proposed by Grimme [2IH21] • This method relies on corrections added 
to the DFT total energy and forces, based on a damped atom-pairwise potential CqR~^ {Cq 
represents the dispersion coefficient for a given atom pair and R is the distance between 
the atoms). The studied system is modelled using supercell, which has an (11 x 11) lateral 
periodicity and contains one layer of (12 x 12) graphene on four-layer slab of metal atoms. 
Metallic slab replicas are separated by ca. 18 A in the surface normal direction, leading to 
an effective vacuum region of about 15 A. To avoid interactions between periodic images of 
the slab, a dipole correction is applied [25]. Due to the large lateral periodicity the surface 
Brillouin zone is sampled with a single fc-point at P for structure optimisation (the positions 
of the carbon atoms as well as those of the top three layers of Rh are optimised) and set to 
3x3x1 mesh in the total energy calculations. The Rh-Rh spacings in the bottom layer 
are hxed at the optimised bulk value. The STM images are calculated using the Tersoff- 
Hamann formalism [2B], which states that the tunnelling current in an STM experiment is 
proportional to the local density of states (LDOS) integrated from the Fermi level to the 
bias. In the visualisation software Hive 1271 it is implemented in its most basic formulation, 
approximating the STM tip by an inhnitely small point source. The integrated LDOS is 
calculated as p(r, e) oc p(r, e')de' with Ep the Fermi energy. An STM in constant current 
£ mode follows a surface of constant current, which translates into a surface of constant 
integrated [p(x, y, z,e) = C with C a real constant]. For each G, this construction returns a 
height X as a function of the position {x,y). This height map is then mapped linearly onto 
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a corresponding colour scale. In AFM images simulations, the tip-sample force is expressed 
as a function of the potential Vts (r) on the tip due to the sample: Fts oc — V[|Vhts(i’)P] 
and transferred to the frequency shift of the oscillating sensor via formalism described in 
Ref. [2H]. 

The graphene/Rh(lll) system was prepared in ultra-high vacuum station for STM/AFM 
studies via cracking of propylene gas (CaHg) according to the recipe described in details in 
Refs. [TTl [121 l29]- Prior to the graphene preparation the Rh(lll) was cleaned via cycles 
of Ar+-sputtering and high-temperature annealing until the clean atomically-resolved im¬ 
ages were acquired in STM. Further, the quality and homogeneity of the graphene/Rh(lll) 
samples were verihed by means of low-energy electron diffraction (FEED) and STM. The 
STM/AFM images were collected with Aarhus SPM 150 equipped with KolibriSensor"*"^ 
from SPECS [301 El] with Nanonis Control system. In all measurements the sharp W- 
tip was used which was cleaned in situ via Ar’''-sputtering. In presented STM images the 
tunnelling bias voltage, Ut, is referenced to the sample and the tunnelling current. It, is 
collected by the tip, which is virtually grounded. During the AFM measurements the sensor 
was oscillating with the resonance frequency of /o = 1001541 Hz and the quality factor of 
Q = 32323, and the frequency shift was used as an input signal in a feedback loop for the 
topography measurements. The oscillation amplitude was set to A = 300 pm. The system 
base pressure was better than 8 x 10“^^ mbar during all experiments. All measurements 
were performed at room temperature. 


III. RESULTS AND DISCUSSION 

The graphene/Rh(lll) can be considered as a representative example of the lattice- 
mismatched graphene-metal interfaces [Fig. [^a)]. In this structure a graphene layer with 
(12x 12) periodicity is arranged on (11x11) Rh(lll) slab. Several high-symmetry adsorption 
sites for carbon atoms in a graphene layer can be identihed in this structure depending 
on the stacking of underlying Rh: ATOP [A] carbon atoms are placed above Rh(S-l) 
and Rh(S-2) atoms], HOP [H; carbon atoms are placed above Rh(S) and Rh(S-2) atoms], 
FCC [F] carbon atoms are placed above Rh(S) and Rh(S-l) atoms], and BRIDGE [B] 
Rh(S) atoms bridge the carbon atoms]. They are marked in Fig. [^a) as circles, triangles, 
squares, and stars, respectively. Obviously that the local symmetry in the graphene/Rh(lll) 
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FIG. 1: (a) Ball model of the graphene/Rh(lll) system. Circles, triangles, rectangles, and starts 
denote the high-symmetry positions for carbon atoms on Rh(lll). (b,c) Top and side views, 
respectively, of the difference electron density maps, An(r) = ngr/Rhif) ~ ’nRhir) ~ ploted 

in units of e/k? calculated for graphene/Rh(lll). 
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structure will define the strength of interaction between graphene and Rh(lll) and among 
all high-symmetry positions the BRIDGE places are expected to be the most energetically 
favourable for the nucleation of the adsorbed atoms and molecules. 

The model shown in Figj^a) was used in the DFT-D2 structure-optimization procedure 
in order to get an information about the spatially distributed interaction strength in the 
graphene/Rh(lll) system. The results of this step are summarised in Fig. [^b,c) where 
the top and side views of the optimised structure of graphene/Rh(lll) together with the 
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distribution of the difference electron density, An(r), are shown. Carbon atoms in the 
ATOP positions dehne the largest distance between graphene and Rh(lll) of 3.15 A. The 
lowest position in this structure is the BRIDGE one with the distance between graphene 
and Rh(lll) of 2.08 A giving a maximum corrugation of the grapehene layer of 1.07 A. The 
HOP and FCC high-symmetry positions are lying by 0.4 A and 0.8 A, respectively, higher 
than the minima. 

The obtained results point out the importance of the dispersive (van der Waals) forces 
for the description of the strongly-corrugated lattice mismatched graphene/metal interfaces. 
The test calculations performed without inclusion of the van der Waals interactions (stan¬ 
dard DFT-PBE approach) give the larger corrugation of 1.8 A with a similar graphene/Rh 
distance for the BRIDGE positions. At the same time the highest position for carbon atoms 
is lying at 3.90 A. This is due to the fact that standard GGA calculations give a reasonable 
result when describing the strong interaction between graphene and metal, whereas the dis¬ 
persion forces, neglected by the standard procedure, are important for the description of the 
weakly interacting graphene/metal interfaces. The latter fact was shown to be valid when 
describing the weakly adsorbed graphene layer on Ir(lll), where bonding was explained by 
the existence of attractive van der Waals forces between graphene and Ir that modulated 
by weak bonding interactions at the EGG and HGP places and anti-bonding interaction 
around ATOP positions |T5]. 

Figs. Bb. c) show the top and side views of the graphene/Rh(lll) interface with the 
difference electron density map, An(r): red colour - electron charge accumulation and blue 
colour - electron dehciency. One can clearly see that around the strongly bonded regions 
the electron density transfer from metal to graphene is observed: Rh 5sp electrons are more 
mobile compared to 4d and the transfer of these electrons on the graphene-derived tt* states 
dehne the doping level of whole graphene which is in metallic state in this system. This effect 
is also rehected in the shift of the corresponding carbon-atom-projected density of states in 
the electronic structure of graphene/Rh(lll): the graphene-derived states are shifted by 
~ 1.5 eV to the larger binding energies compared to the free-standing graphene (Fig. |^. 
Additionally the sublattice symmetry for two carbon atoms in the graphene layer is broken 
at all places. However, only around the strongly-interacting regions it is accompanied with 
the strong hybridisation of the graphene vr states and Rh Ad states that leads to the opening 
of the energy gap for the tt states around the K point of the Brillouin zone of graphene. This 
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FIG. 2: Carbon atom-projected total density of states (cr and tt) in the valence band for the different 
high-symmetry positions of the graphene/Rh(lll) system. The inset shows the corresponding 
density of states for the pz character only. 



effect of hybridisation leads to the appearing of the so-called interface states which can be 
recognised in the DOS picture as a number of peaks around Ep (Fig. |^. The similar effects 
were also observed for the strongly interacting graphene/Ni(lll) and graphene/Fe/Ni(lll) 
systems [321135] . 

We have studied the temperature dependence of the graphene growth at different tem¬ 
peratures and the same pressure of the propylene gas (p = 2 x 10“^mbar). The results are 
compiled in Fig. The growth of the ordered carbon layer (graphene) undergoes several 
stages where disordered carbon layer [Fig. |^b)] or mixture of several carbon phases (carbidic 
and graphene) are observed [Fig. [^c,d)]. Initially, the low-temperature cracking of propylene 
(T = 600 K) on the clean Rh(lll) surface [Fig. |^a)] leads to the formation of the carbon 
layer without any long-range order [randomly distributed carbon-bubbles in Fig. [^b)] and 
only short-range order of carbon rings is observed [inset of Fig. [^b)]. Increasing of the 
synthesis temperature [T = TOOK and T = BOOK for Fig. c) and (d), respectively] leads 
to the ordering of the carbon layer and two carbon phases can be clearly distinguished [A 
(C) and B (D) on the corresponding STM images]. We assign the hrst phase [A in (c) and 
C in (d)] to the carbidic phase of carbon on Rh(lll) and the second phase [B in (c) and D 
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FIG. 3: (a) STM image (150 x 150 nm^) of the clean Rh(lll) surface. The inset (1.9 x 1.9nm^) 
shows the atomically resolved zoom of (a), (b-d) Carbon layers prepared on Rh(lll) at different 
cracking temperatures: (b) 600K, image size 100 x 100nm^, inset size 8 x 8nm^; (c) TOOK, image 
size 70 X 70 nm^, inset size 2.8 x 2.8 nm^ of region A; (d) 800 K, image size 80 x 80 nm^, inset size 
10 X 10 nm^ of the border region between regions C and D. 



in (d)] to the graphene phase on Rh(lll). Note: the studies of the process of the carbon 
phases on Rh(lll) is out of the scope of the present manuscript. However, we would like to 
emphasise that the assignment of the later phases to the ordered graphene layer on Rh(lll) 
is made on the facts that this phase starts to prevail at higher synthesis temperature and 
the observed STM images of this phase are in very good agreement with the calculated ones 
(see discussion below). 

The further increasing of the cracking temperature to T = 1100 K results in the formation 
of the high quality single domain graphene layer on Rh(lll) [Fig. |^a)]. The corresponding 
FEED image of the obtained sample is shown as an inset. The zoomed STM image of the 
graphene moire structure on Rh(lll) is shown in Fig. |^b) where one can clearly identify 
all high symmetry positions discussed earlier (these regions are marked by the respective 















FIG. 4: (a) Large scale STM image (150 x 75 nm^) of the single domain graphene layer on Rh(lll). 
The inset shows the corresponding LEED image, (b) Atomically resolved STM image (6 x 6nm^) 
of the moire structure of the graphene/Rh(lll) system. Tunnelling conditions: Ut = —0.55 V, 
It = 10 nA. (c) Calculated STM image of graphene/Rh(lll): integration was performed in the 
energy range E — Ep = —0.5eV and the distance between tip and the sample was set to 2 A. 



capital letters). The obtained results are well reproduced in the theoretically calculated 
STM image of graphene/Rh(lll) which is presented in Fig. |^c): the corrugation as well as 
the difference in hight for different high-symmetry places are well conhrmed. These results 
compiled on this hgure are in very good agreement with earlier published data and the 
detailed discussion and interpretation of STM data can be found elsewhere [TOHT^ . 

The variation of the charge distribution [Fig. Sb, c)] and the corresponding local density 
of states [Fig. in the graphene moire on Rh(lll) gives a hint on the expectation of the 
most energetically favourable places of this structure for the adsorption and nucleation of the 
deposited atoms and molecules. With this respect the HCP, FCC, and BRIDGE positions 
are the most probable places due to the increased local density of states at Ep for these 
places as well as to the localisation of the negative charge around them. That can lead to the 
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formation of the bonds (physisorption and chemisorption) between graphene and deposited 
atoms. In order to prove this idea we performed the atomic force microscopy/spectroscopy 
experiments on graphene/Rh(lll). Here the oscillating tip can mimic the deposited metallic 
cluster and the interaction force between graphene and tip is reflected in the change of the 
resonance frequency of the sensor. 

The results of these experiments and modelings are compiled in Figs. and The 
resonance frequency shift of the oscillating scanning tip as well as the calculated force and 
the simultaneously measured tunnelling current as a function of the distance between tip 
and sample are shown in the upper row of Fig. The zero-level for the distance is taken 
with respect to the scanning position of the tip when the STM feedback loop is switched on 
to keep tunnelling current constant. One can clearly see that the strength of the interaction 
of tip and sample is different for for ATOP (curve A) and BRIDGE (curve B) positions as 
can be concluded from the discussion placed earlier. At the same time the tunnelling current 
measured in the same distance-dependent experiment shows the exponential dependence as 
expected. 

Figs. ^ a-d) show the constant-frequency AFM images of graphene/Rh(lll) acquired at 
set points marked by the dashed lines and the corresponding letters in the A/ plots from the 
upper row. The respective high-symmetry places for carbon atoms of the graphene/Rh(lll) 
structure are marked by the capital letters. One can clearly see that increasing of the 
frequency shift set-point during scanning leads to the increasing of the imaging contrast. 
Initially, the contrast between hills and valleys can only be distinguished. The further in¬ 
creasing of the frequency shift helps to resolve not only moire structure of graphene/Rh(lll), 
but also to get the atomic contrast in AFM images. This can be explained as due to the 
large contribution of the chemical forces (attractive contribution) when the frequency shift is 
increased. It is interesting to note that the constant frequency shift AFM imaging contrast 
for the graphene/Rh(lll) is exchanged when one going from Figs, [^c) to (d). This effect 
can be understood on the basis of the A/ curves (upper row of Figs. which intersect 
at A/ —lA5Hz indicating different z set points when measuring at the frequency shifts 
corresponding to Figs. |^c) and (d). 

We also perform simulations of constant frequency shift and constant hight AFM images 
of graphene/Rh(lll) (see discussion of the simulation details in the experimental part and 
in Ref. [2S]). Generally, the model analysis of the corrugated graphene-based systems was 
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FIG. 5: (Upper row, from left to right) Frequency shift plots (A/), corresponding interaction force 
(F), and the tunnelling current between tip and sample as functions of the relative distance with 
respect to the scanning position {Ut = —840meV, It = 0.48 nA) for ATOP (A) and BRIDGE 
{B) positions of the graphene/Rh(lll) system, (a-d) Constant frequency shift AFM images of 
graphene/Rh(lll) measured at the frequency shifts of the sensor marked by the dashed lines and 
the corresponding letters in the figure for A/ in the upper row. 



performed in Ref. [36] where the Lennard-Jones model potential between tip and sample was 
taken instead of the real interaction potential. Two contributions in the AFM imaging were 
considered: (i) due to the corrugation of the graphene-based system [in our case corrugation 
of graphene on Rh(lll), i. e. geometrical effect] that shifts two interaction curves in z 
direction and (ii) different strength of interaction between tip and sample for two different 
places as an effect of different interaction of graphene with a substrate at different places 


11 

























[in our case that can be connected with the different local density of states in graphene on 
Rh(lll), i. e. electronic effect] that shifts two interaction curves in vertical direction. 

The results of these simulations are presented in Fig. where they are compared with 
the atomically resolved constant frequency shift AFM image of graphene/Rh(lll). One 
can clearly see very good agreement between two AFM images, experimental and theoret¬ 
ical [Fig. |^a,b)] where the correct corrugation as well as details of AFM images are well 
reproduced. For the constant hight AFM images it is interesting to see that the imaging 
contrast for this structure is fully reversed with respect to constant frequency shift AFM 
image. This can be also understood on the basis of consideration of the attractive parts of 
the A/ curves for different places of the graphene/Rh(lll) structure (Figs.j^: The ATOP 
(A) positions will give more negative frequency shift compared to the BRIDGE (B) posi¬ 
tions (in terms of the absolute values the A/ value is larger for A position). In case of the 
repulsive regime one can expect the opposite situation and the imaging contrast is reversed 
with respect to the attractive imaging regime. Unfortunately, the used approach can not 
give an information about absolute values of the frequency shift and the corrugation of the 
system (see Ref. [2H])- Moreover, the information about repulsive regime in AFM can not 
also be obtained from these simulations. However, if the model repulsive potential is added 
then the full simulation of the system can be performed. 


IV. CONCLUSION 

The geometry and electronic structure of the graphene/Rh(lll) system were studied via 
combination of DFT, STM, and NC-AFM methods. We perform simulation of STM and 
AFM images on the basis of the optimised structure of graphene/Rh(lll) and the charge 
distribution in this system and found that in both methods the imaging contrast is mainly 
defined by the structural corrugation of the system. The simulation of the imaging contrast 
in AFM was performed for the constant frequency shift and the constant hight methods and 
we found that it is inverted with respect to each other when forces between tip and sample 
are attractive. The obtained results shed light on the interactions in the graphene/Rh(lll) 
system and help us to understand the observed effects in STM and NC-AFM imaging of the 
corrugated graphene-based systems. 
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FIG. 6: (a) Experimental constant frequency shift (CFS) AFM image of graphene/Rh(lll). (b) 
and (c) Calculated CFS and constant hight (CH) AFM images of graphene/Rh(lll). 
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